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Abstract. Motivated by recent hints for sterile neutrinos from the reactor anomaly, we study 
active-sterile conversions in a three-flavor scenario (2 active -|- 1 sterile families) for three dif- 
ferent representative times during the neutrino-cooling evolution of the proto-neutron star 
born in an electron-capture supernova. In our "early model" (0.5 s post bounce), the Ve-^s 
MSW effect driven by Am^ = 2.35 eV^ is dominated by ordinary matter and leads to a com- 
plete Ve-Vg swap with little or no trace of collective flavor oscillations. In our "intermediate" 
(2.9 s p.b.) and "late models" (6.5 s p.b.), neutrinos themselves significantly modify the 
Ve-^s matter effect, and, in particular in the late model, vv refraction strongly reduces the 
matter effect, largely suppressing the overall Ve-^s MSW conversion. This phenomenon has 
not been reported in previous studies of active-sterile supernova neutrino oscillations. We 
always include the feedback effect on the electron fraction Yf> due to neutrino oscillations. 
In all examples, is reduced and therefore the presence of sterile neutrinos can affect the 
conditions for heavy-element formation in the supernova ejecta, even if probably not enabling 
the r-process in the investigated outflows of an electron-capture supernova. The impact of 
neutrino-neutrino refraction is strong but complicated, leaving open the possibility that with 
a more complete treatment, or for other supernova models, active-sterile neutrino oscillations 
could generate conditions suitable for the r-process. 
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1 Introduction 

Sterile neutrinos are hypothetical gauge-singlet fermions that could mix with one or more of 
the active states and thus show up in active-sterile flavor oscillations. Low-mass sterile neu- 
trinos have been invoked to explain the excess z?e events in the LSND experiment [1-3] as well 
as the MiniBooNE excess events in both neutrino and antineutrino channels. Interpreted in 
terms of flavor oscillations, the MiniBooNE data require CP violation and thus no less than 
two sterile families [4-6]. However, these models show tension with other oscillation data 
and may require additional ingredients such as non-standard interactions [7]. Moreover, part 
of the parameter space has been excluded by IceCube data [8]. On the other hand, the cos- 
mic microwave background anisotropies [9-13] and big-bang nucleosynthesis [14, 15] suggest 
cosmic excess radiation compatible with one family of sub-eV sterile neutrinos. However, for 
eV-mass sterile neutrinos to be cosmologically viable, additional ingredients are required [16] . 

Our study is motivated by the most recent indication for the possible existence of eV- 
mass sterile neutrinos coming from a new analysis of reactor i/g spectra and their distance 
and energy variation [17-19]. The data suggest a Ve-i^s mixing of sin^ 26 ~ 0.14 and a mass 
splitting of Am^ > 1.5 eV^. In the supernova (SN) context, these parameters imply that 
the Ue flux would undergo MSW conversions to I's closer to the SN core than any other 
oscillation effects. (We assume that, because of cosmological neutrino mass limits, the sterile 
state is heavier than the active ones so that the MSW effect would occur between and 
I's and not between i/g and i's-) Even then, however, the conversion probably would not be 
close enough to the SN core to affect shock reheating during the accretion phase. Moreover, 
removing the flux would stabilize the remaining flux of active neutrinos against collective 
flavor conversions which anyway are likely irrelevant for shock reheating [20-23]. At a larger 
radius, the lost fg flux would be partly replenished by active-active MSW oscillations, in detail 
depending on the mixing parameters among active neutrinos. So while the i^g flux arriving 
at Earth from the next nearby SN would be significantly modified, observational signatures 
would probably require a large fg detector, in contrast to the existing large detectors that 
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primarily measure the arriving Vf, flux by inverse beta decay. Still, possible observational 
signatures may deserve a dedicated study. 

We here focus on a different aspect of Vf.-Vs oscillations that could have an interesting 
impact during the neutrino-cooling phase of the proto- neutron star (PNS). The neutrino- 
driven matter outflow contributes to SN nucleosynthesis and, in particular, it is a candidate 
site for the formation of elements beyond iron by the rapid neutron-capture process (for a 
review, see [24] and references therein). This "r-process" requires a neutron- rich environ- 
ment, i.e. an electron fraction per baryon Yf. < 0.5, sufficiently large entropy to favor a high 
ratio of free neutrons to "seed" nuclei (the latter are usually iron-group nuclei reassembled 
from free nucleons and acting as starting points of neutron captures), and sufficiently fast 
timescales to lower the efficacy of converting alpha particles to heavier nuclei. In standard 
SN simulations, these conditions have remained elusive. The idea that removing the flux 
by active-sterile oscillations can favor a neutron-rich outflow environment is not new [25- 
31]. However, the used mass differences were larger and the possible impact of collective 
active-active oscillations [32] was not taken into account. On the other hand, several recent 
studies have considered the role of collective flavor oscillations (without sterile neutrinos) on 
nucleosynthesis processes like the r-process and the i/p-process in SN outflows [33, 34]. 

In our study we explore the impact of Vf^-Vg oscillations on the electron fraction Y^, based 
on a self-consistent SN model and the corresponding flavor-dependent neutrino fluxes from 
a spherically symmetric (ID) simulation of an exploding electron-capture SN, which leaves 
behind a neutrino-cooling PNS [35]. To this end we begin in section 2 with a description 
of our electron capture SN reference model, we define our notation and fix the neutrino 
mixing parameters. In section 3 we describe the evolution in SN outfiows. In section 4 we 
present our results for three representative times after core bounce (t = 0.5, 2.9 and 6.5 s). 
Conclusions and perspectives are presented in section 5. 



2 Input for neutrino flavor evolution in electron-capture supernovae 

Electron-capture supernovae, originating from low-mass progenitors (8-10 Mq), might repre- 
sent up to about 30% of all core-collapse supernovae [36, 37]. We use long-term simulations 
of a representative progenitor with mass 8.8 [35], performed with the equation of state 
of Shen et al. [38]. For the present study we chose Model Sf 21 (see reference [35] for further 
details; the number 21 denotes a recomputation of the published model with 21 energy bins 
in the neutrino transport instead of the standard 17 bins). In the chosen model, the accretion 
phase ends already at ~ 0.2 s post bounce when neutrino heating reverses the infall to an 
explosion, and the subsequent deleptonization and cooling of the PNS take ~ 10 s. In this 
section, we discuss our reference model for an electron-capture SN, the neutrino fiuxes, and 
the fiavor evolution equations. 

2.1 Reference neutrino signal from electron-capture SNe 

At a radius r, the unoscillated spectral number fluxes for flavor (/3 = e, e, x with x = fi 
or t) are 
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Table 1. Parameters for the three considered, representative post-bounce times (in seconds) of our SN 
explosion modeL Neutrino-sphere radius Ri, in km (assumed equal for all flavors), flavor-dependent 
luminosities L^^^ (in 10^^ erg s~^), average energies {E^^^) (in MeV), and the spectral shape factor ap 
arc listed. 



where L^^ is the luminosity for flavor i/^, {Ey^) the mean energy, and fupiE) a quasi-thermal 
spectrum. We describe it schematically in the form [39] 

/.,(£;) = e,(^)"%-(-^+^)^/^^^^^ (2.2) 

The parameter ap is defined by {E^^^) / {Ey^^)"^ = {2 + a(s)/{l + ap) and is a normalization 
factor such that J dE fy^ (E) = 1 . We choose three representative times during the PNS 
cooling: t = 0.5, 2.9 and 6.5 s after core bounce. In table 1, the neutrino-sphere radius, the 
luminosity Ly^, the average energies {Ey^), and the factor a/j are reported. Concerning the 
neutrino emission geometry, we adopt the spherically-symmetric bulb model [40], with the 
neutrino-sphere radii assumed equal for all active flavors. 

2.2 Neutrino mixing parameters and flavor evolution equations 

Our work is motivated by the reactor anti- neutrino anomaly that requires, if interpreted 
in terms of sterile neutrinos Ug, sizable i'e~'^s mixing. For simplicity we will ignore possible 
mixings of Ug with other active flavors. The mass difference would have to be in the eV range, 
so cosmological hot dark matter limits imply that the sterile state would have to be heavier 
than the active flavors. Besides fe~i^s oscillations, we also include active- active oscillations 
driven by the atmospheric mass difference between fg and Ux and the mixing angle 0i3. 
The and I'r fluxes in a SN are very similar and these two flavors play symmetric roles. 
Therefore, it is useful to deflne a linear combination that is essentially identical with the 
mass eigenstate and mixes with I'g by means of the small G13 mixing angle, and another 
combination that mixes through the solar angle G12 and is separated with the solar mass 
difference Snisoi. Oscillations driven by these parameters tend to take place at a larger radius 
and likely do not affect SN nucleosynthesis. Overall, therefore, we study a 3-flavor problem 
consisting of u^, Ux and Vg with the mass splittings 

<5^atm = -2 X 10-3 eV2 , (2.3) 
5ml = 2.35 eV^ , (2.4) 

where the latter is representative for the reactor-inspired values [17]. Note that we assume 
normal hierarchy for the sterile mass-squared difference 5ml > and inverted hierarchy for 
the atmospheric difference, 5ml^^ < for the atmospheric sector. The associated "high" (H) 
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and "sterile" (S) vacuum oscillation frequencies are then 

, , _ ^rnl^ra _ 5.07 ^ 

Jm^ 5.96 X 10^ 1 , , 

whereas the usual "low" (L) frequency corresponds to the solar mass difference. For the 
active-sterile mixing we use 

sin^ 2614 = 0.165 . (2.7) 
We assume a small mixing between the active flavors, 

sin^ 913 = 10"^ . (2.8) 

In this case the MSW effect driven by this mixing angle is non-adiabatic, i.e. we only focus 
on active-sterile MSW oscillations and collective active-active oscillations. Recent hints for a 
not-very-small value for G13 [41] would imply that we also need to include flavor conversion 
by the active-active MSW effect, but this has little impact on our results. 

We treat neutrino oscillations in terms of the usual matrices of neutrino densities pE 
for each neutrino mode with energy E where diagonal elements are neutrino densities, off- 
diagonal elements encode phase information caused by flavor oscillations. Moreover, we 
work in the single-angle approximation where it is assumed that all neutrinos feel the same 
average neutrino-neutrino refractive effect. The radial flavor variation of the quasi-stationary 
neutrino flux is given by the "Schrodinger equation" 

idrPE = [^e-,Pe\ and idrpE = [^e,Pe\ ■, (2.9) 

where an overbar refers to antineutrinos and sans-serif letters denote 3x3 matrices in flavor 
space consisting of i/g, fx and Ug. The initial conditions are pE = diag(n^^, nj^^, 0) and 
pE = diag(np^, ?ip^ , 0). The Hamiltonian matrix contains vacuum, matter, and neutrino- 
neutrino terms 

He = HT + ^e + ^e- (2.10) 

In the flavor basis, the vacuum term is a function of the mixing angles and the mass-squared 
differences 

Hf= = Udiag(-^,+^|i,^s)ut , (2.11) 

where U is the unitary mixing matrix transforming between the mass and the interaction 
basis. The matter term includes both charged-current (CC) and neutral-current (NC) con- 
tributions and it is in the flavor basis spanned by (fe, fs) 

H- = V2G^ diag(iVe - ^, -^,0) , (2.12) 

where N^, is the net electron number density (electrons minus positrons), and Nn the neutron 
density. 

In all neutral media, Y(. = Yp and Yn = 1 — Ye, where Yj is the number density of particle 
species j relative to baryons. The local electron fraction is 
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Inserting the previous expression for in equation (2.12), the matter Hamihonian becomes 



H- = V2GpN, diag - I, \y, - i , 



(2.14) 



where Ni, is the baryon density. The matter potential can be positive or negative. For 
1^ > 1/3 it is z^e that can undergo an active-sterile Mikheev-Smirnov-Wolfenstein (MSW) 
resonance, whereas for 1^ < 1/3 it is i^e [42]. 

The corresponding 3x3 matrix caused by neutrino- neutrino interactions vanishes for all 
elements involving sterile neutrinos [43], i.e. Hg^ = H^^ = H^^ = 0, and only the 2x2 block 
involving the active flavors is non-zero. In particular, the only non-vanishing off-diagonal 
element of the 3x3 matrix is Hg^. 

In summary, the matter plus neutrino-neutrino part of the Hamiltonian has the diagonal 
elements 



H 



ra+uu 



\/2Gp 



= V2G, 



(2.15) 
(2.16) 



whereas initially the off-diagonal elements vanish. These expressions represent the energy 
shift of i^e or relative to Vg caused by matter and neutrino refraction. 



3 Electron fraction evolution 



The material in a fluid element moving away from the SN core will experience three stages 
of nuclear evolution. Near the surface of the neutron star, typically the material is very hot 
and essentially all of the baryons are in the form of free nucleons. As the material flows 
away from the neutron star, it cools. When the temperature T < 1 MeV, a particles begin 
to assemble. As the fluid flows farther out and cools further, heavier nuclei begin to form. 
Around half of the nuclei with masses A > 100 are supposed to be created by the r-process, 
requiring neutron-rich conditions. In this section, we introduce the Yg evolution equation to 
study whether the impact of sterile neutrinos can help to produce such an environment. 

Having in mind the overall evolution of abundances with radius and time, namely that 
close to the neutrino sphere only free nucleons exist, then alpha-particles begin to form and 
afterwards (some) heavy nuclei, the electron abundance introduced in equation (2.13) can be 
expressed as 

h 

Here Xp (Xa) is the mass fraction of free protons (alpha particles) and Zf^ and Afi are the 
charge and mass number of nuclear species h. The summation runs over all nuclear species 
h heavier than a particles. However, at the conditions common to neutrino-heated outflows 
(in particular in the region where neutrino interactions have the biggest impact on Yg), free 
nucleons and alpha particles typically account for most of the baryons. 

The CC weak interactions alter the electron fraction by converting neutrons to protons 
and vice versa. The electron abundance Yg in neutrino-heated material flowing away from 
the neutron star is set by a competition between the rates of the following neutrino and 
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antineutrino capture reactions on free nucleons, assuming that the reactions of neutrinos on 
nuclei are neghgible, 



Ve + n ^ p + e , (3-2) 
Ve + p ^ n + , (3-3) 

and by the associated reverse processes. Because of slow time variations of the outflow 
conditions during the PNS cooling phase, a near steady-state situation applies and the Yf> 
rate-of-change within an outflowing mass element may be written as [44] 

(iV dY 

-dT = ""^'^^ ^ ^^-"^ + ^'^^^^ " ^^'^ + ^'-^^^^ ' ^^-^^ 

where v{r) is the velocity of the outflowing mass element and and Yp are the abundances 
of free nucleons. The forward rates of the neutrino capture processes of equations (3.2) and 
(3.3) are [44] 

The rates for the reverse processes (electron and positron capture rates on free nucleons) are 
approximately [44] 

A,- 1.578 X 10- s- (^) ' e(--— )/-^ (l + '-^^^ + OmMeV^^ 
V 1.578 X 10-2 g-i f g{-o.5ii-p.)/r. 

1.16 MeV 0.601 MeV^ 0.178 Me 0.035 MeV^\ 
1 H j7 \ \ \ j^i J > (3-8) 

where fie is the relativistic electron chemical potential (in MeV). 

The electron temperature profile is extracted from the hydrodynamical simulation for 
Model Sf 21 of [35] and is shown in figure 1 for the considered outflow trajectories at different 
times after core bounce. The electron chemical potential has been computed by inverting the 
equation [45] 

Ye = ^T^^iv' + n'), (3.9) 

where rj = fJ,e/T is the electron degeneracy parameter. 

Though the details are more complex, Yg at small radii is mainly determined by the e~ 
and e"*" capture rates whereas at larger radii the reverse neutrino-capture reactions dominate. 
Note that in equation (3.4) the nucleons involved in the /3-reactions are free. Including the 
corrections due to nucleons bound in a particles, the free proton and neutron abundances 
are 

Yl = Ye-^ and y^f = 1 - - ^ , (3.10) 

where Xa is the mass fraction of a particles. Assuming that nickel-56 is the most abundant 
among the heavy nuclei, we have tested that the correction to Y^ due to nuclei heavier than a 
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Figure 1. Electron temperature as function of radius from the hydrodynamical simulation of Model 
Sf 21 [35] for t = 0.5, 2.9 and 6.5 s after bounce. 
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Table 2. Initial electron abundance at the neutrino sphere i?^ (in km) for our three considered 
post-bounce times (in s). 



particles is negligible. As the electron temperature, the mass fractions of a particles (and of 
heavy nuclei for the test) are taken from the hydrodynamical simulation of Model Sf 21 of [35] . 
Because the dominant abundances are free neutrons and protons and the differences of the 
Ye evolution with and without neutrino oscillations are relatively small, this is a reasonable 
approximation, and a consistent treatment of the composition evolution would yield only 
minor corrections. 

4 Results 

In order to study the impact of sterile neutrinos on Ye and on the neutrino fluxes, we discretize 
the coupled evolution equations (2.9) in the energy range 1-60 MeV and solve them by 
numerical integration together with equation (3.4). The initial conditions for Ye from the 
hydrodynamical simulation are reported in table 2. The un-oscillated neutrino fluxes are 
fixed by the SN model described in section 2 according to the data given in table 1 for 
t = 0.5, 2.9 and 6.5 s. We consider two possible scenarios: one "with neutrino oscillations" 
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Figure 2. Refractive energy difference T4s between and for the snapshot t = 0.5 s. The 
horizontal line marks ±ws for a typical neutrino energy of 15 MeV. 

where we include the dynamical feedback on Y^. due to neutrino oscillations and the other 
one "without neutrino oscillations." 

4.1 Early cooling phase (t = 0.5 s) 

Figure 2 shows the radial profile of the refractive energy difference V^s = H^^'^'^ — H^'^''"^ = 
H'^^'^ between and Vg as given in equation (2.15). This result already includes a self- 
consistent solution for Yf, caused by neutrino oscillations. The horizontal line marks ws for 
a representative neutrino energy of 15 MeV. An MSW resonance arises for both neutrinos 
and antineutrinos close to the neutrino-sphere. However, here the effective potential varies 
extremely fast because of the sudden rise of Y^, preventing any significant flavor conversion 
(extremely non-adiabatic condition). At larger radii, a second resonance occurs for r < 
300 km, but only for neutrinos, causing adiabatic flavor conversion. Therefore, we expect a 
larger flavor modification in the neutrino sector than in the antineutrino one. 

Figure 3 shows the radial profile of the i'e~^x refractive energy difference caused by 
matter alone 



including the self-consistent solution for 1^. In the active-active sector, collective neutrino 
oscillations will occur. Their relevance is not correctly measured by the refractive energy 
difference V^^ = ^/2G■p {N^^ — Np^ — N^^ + Np^) because collective effects are important 
even if this quantity vanishes due to the off-diagonal refractive index. One way to express 
the possible relevance of collective oscillations is in terms of the quantity 



= V2GFYeNb 




(4.2) 




r (km) 

Figure 3. Refractive v^-Vx energy difference caused by matter, and the estimate ji for the 
neutrino-neutrino interaction energy for the t = 0.5 s model. The horizontal line marks cjh for a 
typical neutrino energy of 15 MeV. 

although somewhat different definitions of ^ have been used in the literature. The potential 
/i is plotted in figure 3 as a function of the radius in order to give a comparison of the active 
neutrino abundance with respect to the electron one. For sake of simplicity, Ny^ are fixed 
to their initial values since the correction due to neutrino density variations does not change 
the ratio between and V^. The horizontal line marks cjh for the neutrino energy 15 MeV. 
No MSW resonance in the atmospheric sector is expected and the matter potential is always 
larger than /i. 

In figure 4 we show the spectra for neutrinos (left) and antineutrino (right). The top 
panels show the primary spectra at the neutrino-sphere (no oscillations). Below we show 
the oscillated spectra when matter refraction is included, causing an almost complete MSW 
swap between fe and Vg, but hardly any conversion between and Vg- Finally in the 
bottom panel we show the result after including neutrino-neutrino interactions in the single- 
angle approximation, causing no further modification. After the i'e~^s MSW conversion, the 
e-x difference spectrum is very asymmetric between neutrinos and antineutrinos, essentially 
suppressing collective conversions. Moreover, Fy^{E) < Fy^{E) and the same is true for 
P. The evolution is therefore inhibited because the system is close to a stable equilibrium 
point [46, 47]. Given the initial conditions produced by the active-sterile MSW effect for the 
subsequent collective oscillations, we expect that any multi-angle treatment might produce 
only a smearing of the spectral features without changing the hierarchy among the fluxes of 
different flavors [48, 49]. In fact the large asymmetry between z/e and i^x fluxes (and the same 
for !>) is responsible for inhibiting any possible flavor decoherence effect [50]. 

Figure 5 shows the effect of oscillations on the Ye profile. The MSW swap between 
and I's almost completely removes the original Ve flux and pushes the matter outflow to 
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Figure 4. Spectra for neutrinos (left) and antineutrinos (right) at r = 1000 km in arbitrary units 
(a.u.) for the 0.5 s model. Top: No oscillations. Middle: Oscillated spectra, including only the matter 
effect. Bottom: v-v interactions are also included, but cause no visible difference. 

a more neutron-rich environment, although not far enough to establish obviously favorable 
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Figure 5. Electron abundance for the t = 0.5 s model, with and without oscillations, 
conditions for an r-process. 

4.2 Intermediate cooling phase (t = 2.9 s) 

We next turn to a snapshot during the intermediate coohng phase at t = 2.9 s post bounce. 
In figures 6-10 we show the analogous information as previously for the early cooling phase. 
There are several new effects. One is that the refractive difference between and i^s quickly 
drops in an almost step-like feature, caused by i^g-z^s MSW conversions (figure 6). Another is 
that the t'e-^'x refractive energy difference caused by matter is now much smaller, allowing for 
an MSW effect between the two active flavors in the neutrino sector for the chosen hierarchy 
(figure 7). The neutrino background is responsible for increasing the flux with respect to 
the case with only matter and for averaging out the i>x and i>e fluxes. 

In flgure 10 the electron abundance is plotted as a function of the radius. The Yg proflle 
is lowered compared to the no-oscillation case. In particular, when I'l/ refraction is included, 
the asymptotic Ye value, due to the increase of Ue flux, is more signiflcantly shifted below 0.5 
than at the earlier time of 0.5 s after bounce. 

For comparison, in figure 9 the neutrino and antineutrino fiuxes are shown for larger 
sin^ Oi3. In this case the MSW i'e~i^x conversion is adiabatic and as a result, the matter-only 
oscillations not only lead to an almost complete swap of i^g to z^g, but in addition a partial z/g- 
i/^ swap. Including neutrino-neutrino refraction, since collective effects and MSW resonances 
are occurring all in the same spatial range, a spectral swap with two splits emerges, one in 
the neutrino sector at about 26 MeV, the other in the anti-neutrino sector at about 6 MeV. 
However, the increase of the fiux at high energies counter-balances the increase of the 
Ue fiux in such a way that the resultant Y^, does slightly change but without modifying our 
conclusions. 
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Figure 6. Refractive energy difference Ves between and Ug for the snapshot t = 2.9 s. The 
horizontal line marks ±ws for a typical neutrino energy of 15 MeV. 
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Figure 7. Refractive Vg-i^x energy difference caused by matter, V™, and the estimate fi for the 
neutrino- neutrino interaction energy for the t = 2.9 s model. The horizontal line marks wh for a 
typical neutrino energy of 15 McV. 




Figure 8. Energy spectra for t — 2.9 s as in figure 4. 
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Figure 10. Electron abundance for the t = 2.9 s model for different oscillation cases as indicated. 
4.3 Late cooling phase (t = 6.5 s) 

Finally, we consider the late cooling phase for a snapshot at t = 6.5 s post bounce. In 
figures 11-14 we show the analogous information as in the previous cases. The active-sterile 
MSW effect once more leads to an almost complete fg-i^s swap, if we include only the ordinary 
matter effect. However, once we include refraction, the results change quite dramatically. 
The active-sterile energy difference again drops quickly at a critical radius because neutrinos 
contribute significantly to the matter effects and the MSW conversion shifts the total matter 
effect to zero and then has oscillatory features, apparently causing parametric resonance 
effects in subsequent neutrino oscillations. Flavor conversions differ for each energy mode, 
inducing wiggles in the energy spectra shown on the bottom panels of figure 13. As a result, 
the oscillated spectra are mixed up and, in particular, the Ue and spectra are almost 
coincident. It is interactions that are responsible for repopulating the fg spectrum and, 
consequently increase the Yg value with respect to the case with only matter background as 
visible in figure 14. The asymptotic value drops far below 0.5 when only the matter effect 
is considered, but including I'ly interactions it is pushed back up above 0.5, although in both 
cases Y(, is lowered relative to the no-oscillation case. 
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Figure 11. Refractive energy difference Ves between and i/g for the snapshot t = 6.5 s. The 
horizontal line marks ±ws for a typical neutrino energy of 15 MeV. 
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Figure 12. Refractive lye-i^x energy difference caused by matter, V^, and the estimate jj, for the 
neutrino-neutrino interaction energy for the t = 6.5 s model. The horizontal line marks wh for a 
typical neutrino energy of 15 MeV. 
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Figure 13. Energy spectra for t — 6.5 s as in figure 4. 
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Figure 14. Electron abundance evolution for the t = 6.5 s model for different oscillation cases as 
indicated. 

5 Conclusions 

Motivated by recent hints for the existence of sterile neutrinos, and in particular the antineu- 
trino reactor anomaly, we have studied flavor oscillations within a three- flavor (fe,z^x)i^s) 
scheme in the context of an electron-capture SN model, focussing on the neutrino-cooling 
phase of the nascent neutron star. We have included fe-z^s mixing with parameters suggested 
by the reactor anomaly and active-active mixing representing 1-3 oscillations driven by the 
atmospheric mass difference and a small 0i3 angle. Our main goal was the determination of 
neutrino flux and spectral changes and their impact on the evolution and asymptotic value 
of Ye in the neutrino-driven wind ejecta. Their modification is relevant for nucleosynthe- 
sis in SN outflows. We have studied three snapshots that are representative for the early, 
intermediate, and late cooling stages. 

Even in our simplified neutrino mixing scheme, the results depend sensitively and in 
complicated ways on the detailed matter profile and neutrino fluxes and spectra. In the early 
phase after the onset of the explosion, oscillations are driven almost entirely by the ordinary 
matter effect and lead to a simple i^e-i^s MSW conversion. In the latter cases, neutrinos 
contribute significantly to the refractive energy shifts. One result is that the overall I'e-i^s 
energy difference drops to zero when some of the MSW conversion is complete and in this 
way oscillations feed back on themselves. The switch-off of the matter effect apparently can 
lead to parametric resonance effects and a repopulation of z/g from i^x- 

In all cases, the asymptotic Yg value is lowered compared to the non-oscillation case, 
but it sensitively depends on the cooling phase how large this effect turns out to be. The 
inclusion of active-active collective oscillations and MSW conversions in addition to active- 
sterile mixing strongly modifies the outcome. These general trends are not severely altered 
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by larger values of ©13. 

Accordingly, neutrino conversion to a sterile flavor as well as oscillations and collective 
transformations of active flavors influence the radial variation and time-dependent asymptotic 
value of Y(. in the neutrino-driven wind in complicated ways. The feedback of active-sterile 
oscillations on the refractive effect causes intriguing nonlinear modifications of the naive 
oscillation picture and active-active oscillations can in addition play an important role in 
determining the neutron-to-proton ratio of SN ejecta. 

The chemical composition of the matter outflow can thus be strongly affected by neu- 
trino oscillations. In our model of the neutrino cooling of the proto-neutron star born in 
an electron-capture SN, the corresponding changes do not lead to a large neutron excess. 
Therefore it appears unlikely that in the studied cases viable r-process conditions could be 
produced by flavor oscillations, although the formation of heavy elements must be expected 
to change. 

All numerical studies of neutrino oscillations with neutrino-neutrino refraction use sim- 
plified assumptions, in our case the "single angle approximation," but even more sophisti- 
cated "multi-angle studies" assume axial symmetry around the radial direction and therefore 
are not fully general. Even if we suspect that our conclusions are not affected by the ly-iy 
interaction's angular dependence, a careful analysis on this effect remains certainly to be 
investigated. 

More important may be the calculation of a denser grid of snapshots for more solid 
conclusions on the nucleosynthetic implications and as the basis for detailed studies of element 
formation. Moreover, because the neutrino emission properties and the neutrino-driven wind 
conditions depend sensitively on the mass of the proto-neutron star [51], our results for the 
1.36 Mq (baryonic mass) remnant of an electron-capture SN may be applicable only to SNe 
with similar compact objects. The investigation of a broader range of progenitor models, in 
particular also of iron-core SNe with more massive proto-neutron stars, is therefore desirable 
to identify possible cases where favorable conditions for the r-process may be produced by 
flavor oscillations involving sterile neutrinos. 

If sterile neutrinos with parameters suggested by the presently discussed antineutrino 
reactor anomaly are verified in future experiments, their existence cannot be ignored in 
nucleosynthesis studies of the SN environment. Sterile neutrinos must be expected to have 
important consequences for the possibility of a z^p-process in SN outflows and might even be 
relevant for the question whether SN explosions can be sources of r-process elements. The 
experimental neutrino mixing parameters will be a crucial input information for theoretical 
investigations of these problems. 
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